Universal Long-time Behavior of Nuclear Spin Decays in a Solid 
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Magnetic resonance studies of nuclear spins in solids are exceptionally well suited to probe the 
limits of statistical physics. We report experimental results indicating that isolated macroscopic 
systems of interacting nuclear spins possess the following fundamental property: spin decays that 
start from different initial configurations quickly evolve towards the same long-time behavior. This 
long-time behavior is characterized by the shortest ballistic microscopic timescale of the system and 
therefore falls outside of the validity range for conventional approximations of statistical physics. 
We find that the nuclear free induction decay and different solid echoes in hyperpolarized solid 
xenon all exhibit sinusoidally modulated exponential long-time behavior characterized by identical 
time constants. This universality was previously predicted on the basis of analogy with resonances 
in classical chaotic systems. 
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The relationship between statistical physics and chaos 
is one of the most important and controversial problems 
in theoretical physics. Statistical physics is based on the 
assumption of some kind of randomness on the micro- 
scopic scale, yet the question of whether this randomness 
is at all related to the mathematical concept of chaos 
(well-established for few-body classical systems) is not 
well understood [H, Q. In many-body systems, it is ex- 
tremely difficult to separate the effects of randomness as- 
sociated with true chaos from those associated with aver- 
aging over the macroscopic number of degrees of freedom 
[Hi 0> HI IB]- The situation is further complicated by the 
lack of consensus on the universal definition of chaos in 
quantum systems. In view of these complications, one 
approach is to proceed on the basis of conjectured paral- 
lels between the properties of mathematical chaotic sys- 
tems and real many-body systems. The predicted con- 
sequences of these conjectures can then be tested nu- 
merically or experimentally. One such prediction about 
the universal long-time behavior of transient nuclear spin 
decays in solids has been made recently in Ref. [?]. The 
work presented here tested that prediction by measuring 
the transverse relaxation of 129 Xe nuclei (spin = 1/2) 
in solid xenon over four orders of magnitude using nu- 
clear magnetic resonance (NMR). Such experiments are 
prohibitively challenging for conventional NMR due to 
the weak thermal magnetization achievable in even the 
strongest magnets. We have employed the technique of 
spin-exchange optical pumping [8[ in order to achieve en- 
hanced (hyperpolarized) magnetization required for this 
experiment. 

In nearly perfect agreement with the prediction of 
Ref. 0, our experiments indicate that the long-time be- 
havior of transverse nuclear spin decays in solids has the 
universal functional form 

F(t) = Ae~ 7t cos(wi + <t>), (1) 

where the decay coefficient 7 and the beat frequency u 



are independent of the initially generated transverse spin 
configuration. This long-time behavior sets in after only a 
few times T2, where T2 is the characteristic timescale for 
transverse decay determined by the interaction between 
nuclear spins (see Eq. (J2j) below) and represents the short- 
est ballistic timescale in the system. The values of 1/7 
and l/u> are also on the order of T%. Hence, it cannot 
be that the spins are interacting with a fast-equilibrating 
heat bath, which would justify the exponential character 
of the decay, as for a common damped harmonic oscilla- 
tor. Indeed, at 77 K in an applied magnetic field > 1 T, 
the 129 Xe spins are well isolated from their environment. 
The longitudinal relaxation time T\ « 2.3 h § while 
T2 ~1 ms : therefore, the decay cannot be attributed 
to spin-lattice relaxation. The oscillations in this decay, 
sometimes referred to as Lowe beats [Hj], constitute a 
correlation effect 12, 13| induced by the spin-spin inter- 
action and have nothing to do with the Larmor frequency. 
We verified that the effects of radiation damping and in- 
homogeneities in the external field are also negligible on 
the timescale of T2 [1J] . The observed decay thus repre- 
sents the approach of a closed quantum system to equi- 
librium. 

We used the particular pulse sequence (see Fig. [J), 
known as a solid echo [H, [lj] , which consists of two 90° 
pulses (the first along the y-axis and the second along 
the x-axis in the rotating frame) separated by a delay 
time t. In contrast with the conventional Hahn spin 
echo [l6[ or the magic echo [17], the solid echo is not an 
amplitude-attenuated reproduction of the free induction 
decay (FID) that peaks at time It. Complete refocus- 
ing by solid echoes occurs only for isolated pairs of spins 
IS]. A deviation from complete refocusing is caused by 
higher-order correlations involving more than two spins. 
The solid echo response depends on the spin configura- 
tion just after the second pulse, whereby different values 
of the delay time r imply fundamentally different "after- 
pulse" configurations [l4[ that evolve from the uniformly 
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FIG. 1: The pulse sequence used to generate a solid echo. 
The magnitude of the free-induction decay (FID) including 
Lowe beats is shown schematically after the first pulse with 
the solid-echo response shown after the second pulse. The 
pulses are separated in phase by 90°: the first is along the 
rotating-frame y-axis, the second is along the rotating-frame 
x-axis. Unlike conventional Hahn echoes, the solid echo does 
not generally peak at time 2r[3]. 



polarized uncorrelated spin state at the beginning of the 
FID to highly correlated states induced by spin-spin in- 
teractions during the delay time [l8| . Experimentally, 
these distinct after-pulse configurations are exactly what 
is required in order to clearly demonstrate the evolution 
to a universal long-time behavior. 

On the timescale of our experiments, the system of 
interacting 129 Xe nuclei can be accurately described as 
isolated and governed by the Hamiltonian of the trun- 
cated magnetic dipolar interaction [l9|, [2(| , which in the 
Larmor rotating reference frame has the form 



mi) 



(2) 



where A^ n and Bk n are coupling constants and I % n are 
the spin operators representing the ith projection of the 
nth spin. The Hamiltonian in Eq. ([2]) is appropriate in 
the high-field limit where the Zeeman energy dominates 
dipolar couplings; hence, the shape and duration of the 
FID are independent of the applied field. The character- 
istic decay timescale is on the order of a few inverse 
nearest-neighbor coupling constants. Although the cou- 
pling constants in the Hamiltonian Eq. ([2]) can be very 
accurately determined from first principles, efforts over 
several decades 0, [ll], [H, H3, [U, Hi, H3, [S] to predict the 
entire behavior of FIDs and spin echoes quantitatively 
have met only with limited success; direct calculations 
tend to lose predictive power in the long-time tail of the 
decay, where increasingly higher-order spin correlations 
become important Although methods have not yet 
been developed for the controllable calculation of oj and 
7 in Eq. {T]), Ref. Q predicted the quick onset [ijj of 
the long-time behavior of Eq. with the same values 
of u and 7 for all kinds of transverse decays in the same 
system. This prediction was based not on a conventional 
statistical theory but on a conjecture (l2| that quantum 



spin dynamics generates extreme randomness analogous 
to classical chaos. 

The long-time behavior of Eq. ([1]) for the FID alone 
has been previously observed in NMR experiments on 
19 F in CaF2 [26j]. Analogous observations have also been 
made in numerical simulations of both classical 27J and 
quantum [28| spin lattices. Our experiment adds a new 
insight into this universality by demonstrating that the 
above long-time behavior is common to both FIDs and 
solid echoes in the same spin system. Since solid echoes 
initiated at different delay times r start from distinct 
initial spin configurations, our findings suggest that the 
tails of transient nuclear spin signals are independent of 
initial conditions (apart from the oscillation phase and 
the overall amplitude). 

For our experiments, both isotopically natural (26.4% 
129 Xe, 21.29% 131 Xe) and enriched (86% 129 Xe, 0.13% 
131 Xe) samples of solid polycrystalline xenon contain- 
ing 129 Xe polarized to 5-10% were prepared using spin- 
exchange convection cells [3, HH • FIDs and solid echoes 
were acquired at 77 K in an applied field of 1.5 T ( 129 Xe 
Larmor frequency of 17.6 MHz), well into the high-field 
limit of Eq. Q . The enormous dynamic range of these 
signals required a separate acquisition of the initial and 
long-time decays for each FID and echo using different 
gain settings for the NMR receiver 14 1. 

In Fig. representative decays of the signal mag- 
nitude for the FID and solid echoes with three differ- 
ent delay times t are shown for enriched xenon on a 
semilog plot, with the time axis referenced to « 100 us 
(instrumental dead time) after the end of the first 90° 
pulse, i.e., at the start of the FID. The FID and each 
echo are acquired separately with the sample newly po- 
larized, whereby the run-to-run variation in polarization 
prohibits a direct measurement of their relative ampli- 
tudes. Hence, the data for each echo are shown at the 
proper temporal location, starting w 100 us after the 
corresponding value of r, and each echo is normalized 
to match the FID amplitude at t = r. Fig. shows 
the same four acquisitions time-shifted and amplitude- 
normalized relative to the FID to yield the best overlap 
at long times (« 2.7 ms and later after the start of the 
FID or echo) . The decay coefficient 7 and beat frequency 
to were obtained for each decay from a fit of the long-time 
signal magnitude to the absolute value of Eq. ([1]) ; a repre- 
sentative fit is shown in red. The results are summarized 
in Table HI where each entry represents the average of fits 
for six separate acquisitions of the FID or solid echo. For 
a given isotopic concentration of 129 Xe, the parameters 7 
and u> are the same for both the FID and all solid echoes 
independent of the delay time r. 

In contrast, there is no universal behavior in the ini- 
tial portion of the transverse decays. This is connected 
to the theoretical expectation (discussed above) that the 
longer values of the delay time r allow higher-order spin 
correlations that are not refocused by the solid echo to 
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FIG. 2: Representative acquisitions of FID and solid echoes 
(with three different delay times) recorded for 129 Xe in en- 
riched polycrystalline xenon at 77 K. a. The signal magni- 
tudes are shown on a semilog plot (main) and exhibit char- 
acteristic beats. The actual signals plotted on a linear scale 
(inset) change sign, whereas the main plot has cusps at the 
zero-crossing points. The FID starts at t = ms, and each 
echo starts at its respective delay time r with its initial value 
normalized to the value of the FID at time r. b. The same 
data are shown again on a semilog plot (main) with the echoes 
both time- and amplitude-shifted to illustrate the nearly per- 
fect overlap of the long-time decays. For the 1.9 ms and 2.5 ms 
echoes, the noisiest data farthest out in time have been re- 
moved for clarity. The red line is a representative long-time 
fit to the absolute value of Eq. fT} with the decay coefficient 
7 = 2.04 ms -1 and the beat frequency u = 1.25 rad/ms. The 
distinct differences among the initial portions of the decays 
can be better appreciated in the linear absolute-value plot 
(inset). 



become stronger. As a result, the initial spin configura- 
tions for the FID and various solid echoes are different 
and not trivially related to one another [l4j |. 

In natural xenon, the long-time tails of the FID and 
the solid echo acquired with delay time r = 0.56 ms 
are also nearly identical [3] and can be fit by Eq. (fTJ) 
with parameters 7 and ui given in Table HI The values of 
both of these parameters are smaller than in the enriched 
sample because in the natural sample the 129 Xe spins are 





7 (ms ) 


ijj (rad/ms) 


Enriched FIDs 


1.25 ± 0.05 


2.03 ± 0.04 


T71 'ill r\ f r> 

Enriched echoes, r = 0.56 ms 


1.25 ± 0.05 


2.00 ± 0.03 


Enriched echoes, r = 1.9 ms 


1.22 ±0.04 


2.06 ±0.03 


Enriched echoes, r = 2.5 ms 


1.25 ±0.04 


2.05 ± 0.04 


Natural FIDs 


1.04 ±0.08 


1.53 ±0.08 


Natural echoes, r = 0.56 ms 


1.04 ±0.12 


1.52 ± 0.04 



TABLE I: The decay coefficient 7 and beat frequency w ex- 
tracted from the fit of long-time data by Eq. ((TJ for FID and 
solid echo experiments in both natural and 129 Xe-enriched 
solid xenon. Each entry represents an average of six separate 
experiments with the errors determined from the spread in 
the fit results. The delay time r is the time between the 90° 
pulses in the solid echo pulse sequence. 



more dilute, having been replaced with zero-spin species 
or with 131 Xe, for which the dipolar interaction has dif- 
ferent coupling constants. The intrinsically weaker signal 
meant that the long-time tails of echoes with delay times 
r > 0.56 ms could not be accurately measured. 

The common quantitative long-time character of FIDs 
and spin echoes provides experimental support for the 
notion of eigenmodes of the time evolution operator T(t) 
in isolated many-body quantum systems. This operator 
is defined by the equation £?(i, x) = T(t)g(0,x), where 
g(0, x) is the many-body density matrix at some initial 
time t = 0, and x is the set of variables that describe 
the density matrix. It was conjectured 0, EH that in 
the observable long-time range, the non-equilibrium be- 
havior of the density matrix for any small but macro- 
scopic subsystem of the closed system is controlled by a 
complex- valued eigenmode having the form: 



£>o(x)e 



e5(x) 



(3) 



If this conjecture is valid, then the long-time decay of 
Eq. (fTJ) represents not just the property of one relaxation 
process, such as the FID, but rather an intrinsic property 
of the many-body dynamics of the system, and should 
manifest itself in numerous other relaxation processes, 
such as solid echoes with different delay times r. 

The eigenmodes of the time evolution operator as de- 
fined by Eq. have no direct relation to the eigenvalues 
of the Hamiltonian of the many-body system, but rather 
they are expected to be counterparts of the Pollicott- 
Ruelle resonances 0, [3(J in classical hyperbolic chaotic 
systems. These resonances depend on the rate of proba- 
bility loss from coarser to finer partitions of phase space 
[3,[2|. In many-body quantum systems, there should ex- 
ist an analogous transfer of spectral weight from lower to 
higher order quantum correlations 

Quantum analogs of Pollicott-Ruclle resonances have 
been observed numerically in kicked spin-1/2 chains [3l| . 
the kicked quantum top [32| , Loschmidt echoes [IH , and 
experimentally for the imitation of the single particle 
quantum problem in microwave billiards 3J| . In all these 
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cases, the quantum systems had one or several of the 
following features: (i) very few degrees of freedom; (ii) 
proximity to the classically chaotic limit; (hi) application 
of external time-dependent forces, removing the difficulty 
associated with the discrete frequency spectrum of an iso- 
lated quantum system. In contrast, we deal here with an 
essentially isolated system having a macroscopic number 
of maximally non-classical components (spins 1/2), and 
have no specially introduced precondition for chaos apart 
from the naturally occurring non-integrable interaction 
between spins. 

We note a remarkable fact revealed by Fig. [2^: the 
phases of the long-time oscillations of the 1.9 ms and 
2.5 ms echoes nearly coincide with each other and are 
shifted by it with respect to the FID phase. Indeed, 
one can observe that the zero-crossings (cusps) of the 
FID and the two echoes coincide in the long-time regime. 
Given that these are the absolute-value plots, the above 
coincidences imply that the relative phases of the long- 
time signals are either zero or 7r. These two possibili- 
ties can be discriminated by keeping track of the suc- 
cessive sign changes at the zero-crossings for each signal. 
(The inset of Fig. [2^ shows the sign of the FID and each 
echo.) This may be a fundamental phase relation asso- 
ciated with the fact that the 1.9 ms and 2.5 ms echoes 
start after the FID has begun to approach the asymp- 
totic regime of Eq. ([1]). In contrast, the 0.56 ms echo 
starts well before the FID has reached that regime, and 
its phase has no particular relation to the other three 
signals. 

We have observed a universal long-time behavior of 
129 Xe FIDs and solid echoes in solid xenon. In all cases, 
a sinusoidally modulated exponential decay sets in after 
just a few times T^. This behavior is universal in the 
sense that the two parameters characterizing the long- 
time decay are independent of the NMR pulse/delay se- 
quence, even though each such sequence generates a dif- 
ferent initial spin configuration. These findings reveal a 
fundamental property of nuclear spin dynamics. In addi- 
tion, they also support the idea that the correspondence 
between classical and quantum chaotic properties of real 
many-body systems can be established at the level of 
Pollicott-Ruelle resonances. Further investigations, how- 
ever, are required in order to clarify whether the eigen- 
modes of form (J3J) actually exist in many-spin density 
matrices and, if so, how far this correspondence can be 
taken. 
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SUPPORTING MATERIAL 

(with a separate system of references) 

Detailed Experimental Methods 

The samples of solid polycrystalline hyperpolarized 
129 Xe were prepared using spin-exchange convection cells 
These cells allow the production of large quantities 
of hyperpolarized 129 Xe by keeping most of the xenon in 
the liquid phase. The xenon gas is polarized at 100 °C 
and travels by convection to the column of liquid that is 
kept cold (— 110°C)ina retort that is open to the rest of 
the cell. The liquid 129 Xe is polarized by passive phase 
exchange with the gas. Using these 10 cm 3 cells, we have 
generated a nuclear spin polarization of 5-10% in one liq- 
uid sample of 0.7 mM enriched xenon (86% 129 Xe, 0.13% 
131 Xe ) and two liquid samples of 2.6 mM natural xenon 
(26.4% 129 Xe, 21.29% 131 Xe) inside a 1.5 T horizontal- 
bore superconducting magnet ( 129 Xe Larmor frequency 
of 17.6 MHz). 

To study the NMR signals in the solid state, the liq- 
uid hyperpolarized xenon was subsequently frozen and 
maintained at 77 K while still inside the magnet by forc- 
ing liquid nitrogen directly into the chamber containing 
the cell. At 1.5 T, the polarization is expected to survive 
this phase transition intact 0, and this was verified ex- 
perimentally by monitoring the NMR signal height with 
small-angle pulses during freezing. Since 90° rf excita- 
tion pulses (pulse length w 10 /is) are used to obtain 
the FIDs and solid echoes, the magnetization was es- 
sentially destroyed by each acquisition. However, the 
convection cells allow easy in situ regeneration of the 
129 Xe polarization, and the experiment could thus be re- 
peated many times per day. The NMR flip angle was 
calibrated for each separate measurement by applying 
several identical excitation pulses (8 « 40°) and mea- 
suring the corresponding cos 9 magnetization loss. The 
separate acquisitions of the initial and long-time portions 
of the decays for each FID and echo were joined together 
at one point common to both signals and amplitudc- 
renormalized. (The renormalization is necessary because 
the 129 Xe polarization obtained by spin-exchange opti- 
cal pumping varies somewhat from run to run.) In one 
case (the 0.56 ms delay time in enriched Xe), three sepa- 
rate signals had to be acquired, and these were similarly 
joined together at two locations. 

We have examined the possible effects of radiation 
damping and inhomogeneity in the external field on the 
FID and solid echo line shapes. The observed trans- 
verse decay time Tjf for hyperpolarized liquid xenon (hav- 
ing approximately the same polarization as the solid) is 
« 23 ms in the same NMR probe and location in the 



magnet. (Here, T 2 * refers to the timescale for transverse 
decay from all sources of dephasing combined.) Unlike 
the solid case, the spin-spin interactions in the motion- 
ally narrowed liquid are much weaker, allowing the decay 
to be essentially limited by some combination of radia- 
tion damping and field inhomogeneity. (Since we observe 
the liquid lineshape to be independent of polarization, it 
is likely that field inhomogeneity is actually the limiting 
factor.) Since the liquid decay is much longer than the 
w 1 ms decay time observed for the solid, we conclude 
that the effects of these other sources of dephasing are 
negligible in the solid. As an additional check, we note 
that the decay shapes associated with both the field in- 
homogeneity and the radiation damping are expected to 
be non-exponential; if these factors were not negligible 
in our experiments, then the long-time decay of the beat 
amplitude would exhibit a noticeable departure from an 
exponential shape. 

Theoretical facts about solid echoes 

Solid echoes are produced by the radio frequency pulse 
sequence 90°— time delay r— 90° . Here and below all axes 
are defined in the Larmor rotating reference frame. The 
pulses are assumed to rotate the spin system instanta- 
neously. 

Before the first pulse, the spin system is characterized 
by the equilibrium density matrix 

~k^T J ' ( ' A ' 4 " ) 

where M z — gJ2k ^k 1S t ne total magnetization of nuclei 
with gyromagnetic ratio g in external magnetic field H 
directed along the z-axis, ks is Boltzmann's constant, 
T the temperature (in our case, it is the effective spin 
temperature of the hyperpolarized system), and Cat is 
the normalization constant. The usual energy hierarchy 
ksT 3> JigH 3> h/T-2 is satisfied in our experiment. 
The operator of the first (90°) pulse is 




(A.5) 



This operator transforms If — > P^I^Pi = and, there- 
fore, just after the pulse 

(A.6) 

where M x = gj^k 1 ^- 

Following the first pulse, the system undergoes free 
induction decay (FID) described as follows: 

F FlD (t) - (M x (t)) =Tr{E fc ^e mt p(0) e - mt } 

= Tr{£ fc ^e mt £„^e- mt }, (A . 7) 
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where the second sign = implies that we expanded p(0) 
according to Eq. (|A.6[) , dropped the term containing 1 as 
it makes no contribution to the trace, and then omit- 

C a 2 

ted the prefactor fc ^y • ^ ne Hamiltonian H is given by 
Eq.(2) of the main article. 

The second (90°) pulse is characterized by operator 



P 2 =exp j-ij^j. 



(A. 



The action of this operator is the following: 



Ik 
It 



it 



P]-I V k P2 = 'I 
P}l Z k P2 =lt 



(A.9) 



The density matrix just before the second pulse is 



P-(r) 



p(0) e 



-iHr 



Just after the pulse, it becomes 

p + (r) = P 2 + e ffl V(0)e-'*F, 



(A.10) 



(A.ll) 



The action of P 2 + ...Pi in Eq. ljA.lip can be represented 
as the change of the identity of I k operators according to 
rules (|A.9|) . (Here a = x, y, z.) Indeed, since P2P2 = 1, 



any term of the form P% I k Ii ---In Pi arising in the ex- 
pansion of the right-hand-side of Eq. ljA.lip is equal to 

p 2 Ik p z p 2 T i P2P2-P2P2 !n p 2, and hence each oper- 
ator I? is being "rotated" individually. Using the above 
fact, Eq. (|A.llj) can be rewritten as 



P+ (T) = e m * T p(Q)e 



(A.12) 



where 



H R = £ B kn F k Il + A kn Ipy + B kn P k I*. (A. 13) 



(The subscript "R" stands for "Rotated".) 

The density matrix at time t' after the second pulse is 



psE{T,t) = e P+(T)e 
Finally, the solid echo signal is given by 



(A.14) 



F SE (r,*') = Tr 



x AH t' AHr t 



(A.15) 

As mentioned in the main article, different delay times 
t imply fundamentally different initial spin configura- 
tions, p+{r), for the solid echo response. At small r, 
P+(t) is not much different from p(0) given by Eq. (|A.6[) . 
which describes uncorrelated uniformly polarized spin 
distribution. As r increases, the density matrix p+{r) 
is gradually overtaken by increasingly higher order spin 



correlations induced by the Hamiltonian H R . This can 
be seen from the expansion 

p+(r) = P (0) + ir[H R , J2 k Ik] ~ £ fe Ifi] 

[Wr, [H r , [Wr,E*^]]] + - ■ ( A -!6) 

Each commutation with H R entails an extra power of 
operators I k . The resulting higher order spin correlations 
dominate in p+(r) at r > T2. We further note that in this 
regime, the initial spin configurations for a solid echoes 
with different delay times r are not equivalent to each 
other. 

The correlations induced by H R are non-trivially con- 
nected with those intrinsic for H, the actual Hamiltonian 
of the system. One can achieve a limited simplification 
of the problem by expressing H R as 



H R 



-H + AH, 



(A.17) 



where 



AH 



2B kn I%Il + (A kn +B kn )(qiy + F k P n ) 



(A.18) 

The advantage of this representation is that AH com- 
mutes with Efe I k - If it also commuted with H, then the 
exponent e i (~ w + A ' H ) T would factorize into e~ in T e lAn r , 
and then e lAHT would not influence p+(r), while e~ lHT 
would lead to a time-reversed evolution and imply the 
full recovery of the initial state at t' = r [i.e. p se (t, t) = 
P(0)]- 

As mentioned in the main article, such a situation is, 
indeed, realized for a pair of interacting spins 1/2. In this 

case, [H,AH] = 0, because [I™ I-?, = f° r an y a 

and /?. However, for a lattice of spins 1/2, [H, AH] ^ 0, 
because the commutators like [Jf > I2I3] are n °t equal 
to zero for a ^ (3. 

One can nevertheless conclude that solid echo tends to 
undo two-spin correlations in many-spin systems. There- 
fore, at sufficiently small t, when p+{r) is controlled by 
the first three terms in expansion (|A.16|) . the solid echo 
can be described as an imperfect attempt at time rever- 
sal. (Note that even for small delay time 0.56 ms, the 
maximum of the solid echo in Fig. 2a of the main ar- 
ticle occurs noticeably earlier than at t = 2r.) When 
t <L T2, the rigorous calculation of solid echoes becomes 
intractable, but otherwise it is clear the concept of time 
reversal becomes increasingly inappropriate. 

For further reading, see the original references concern- 
ing solid echoes 0, |i| . 
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Natural Xenon Data 




■ FID (Natural Xe) 

■ 0.56 ms echo (Natural Xe) 
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FIG. SI: Representative acquisitions of FID and the 0.56 ms 
solid echo recorded for 129 Xe in polycrystalline xenon with 
natural abundance of nuclear isotopes at 77 K. Here, the FID 
and echo and normalized and time-shifted in the same manner 
as for enriched xenon in Fig. 2b (main manuscript) to show 
the overlap in the long-time regime. The long-time behavior 
of the FID and solid echo are identical in form, although the 
decay coefficient 7 and beat frequency ui each have values 
smaller than those for enriched xenon; see Table I of main 
manuscript. 



